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Abstract 

The  mechanical  response  of  fuel  cell  proton  exchange  membranes  subjected  to  a  single  hygro-thermal  duty  cycle  in  a  fuel  cell  assembly  is 
investigated  through  numerical  means.  To  this  end,  the  behavior  of  the  membrane  with  temperature  and  humidity  dependent  material  properties  is 
simulated  under  temperature  and  humidity  loading  and  unloading  conditions.  The  stress-evolution  during  a  simplified  operating  cycle  is  determined 
using  finite  element  analysis  for  two  clamping  methods  and  two  alignments  of  the  bipolar  plates.  It  is  shown  that  compressive,  plastic  deformation 
occurs  during  the  hygro-thermal  loading,  resulting  in  tensile  residual  stresses  after  unloading.  These  residual  in-plane  stresses  in  the  membrane 
may  explain  the  occurrence  of  cracks  and  pinholes  in  the  membrane  under  cyclic  loading. 

©  2006  Published  by  Elsevier  B.V. 
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1.  Introduction 

The  durability  of  proton  exchange  membranes  (PEMs)  used 
in  fuel  cells  is  a  major  factor  in  the  operating  lifetime  of  fuel  cell 
systems.  Durability  is  a  complicated  phenomenon,  linked  to  the 
chemical  and  mechanical  interactions  of  the  fuel  cell  compo¬ 
nents  -  i.e.,  electrocatalysts,  membranes,  gas  diffusion  layers, 
and  bipolar  plates  -  under  severe  environmental  conditions,  such 
as  elevated  temperature  and  low  humidity  [1-6].  In  fuel  cell 
systems,  failure  may  occur  in  several  ways  such  as  chemical 
degradation  of  the  ionomer  membrane  or  mechanical  failure  in 
the  PEM  that  results  in  gradual  reduction  of  ionic  conductivity, 
increase  in  the  total  cell  resistance,  and  the  reduction  of  volt¬ 
age  and  loss  of  output  power  [1-3,6].  Mechanical  damage  in 
the  PEM  can  appear  as  through-the-thickness  flaws  or  pinholes 
in  the  membrane,  or  delamination  between  the  polymer  mem¬ 
brane  and  gas  diffusion  layers  [3,5-7].  Although  the  initiation 
and  propagation  of  flaws  and  pinholes  is  not  well  understood,  it 
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has  been  suggested  that  the  stresses  in  the  membrane  are  due  to 
the  hygro-thermal  cycles  even  in  the  absence  of  electric  poten¬ 
tial  or  reactive  gases  [1,2, 4-8].  Thus,  the  purpose  of  the  current 
work  is  to  investigate  the  evolution  of  mechanical  stress  due  to 
the  environmental  conditions  of  a  duty  cycle. 

Work  concerning  testing,  modeling  and  mechanical  charac¬ 
terization  of  perfluorosulfonic  acid  (PFSA)  membranes,  inves¬ 
tigating  the  failure  mechanisms  in  the  membranes  and  how 
hygro-thermal  stresses  influence  in  these  failures,  is  available 
in  the  literature  [1,2,7-12].  However,  failure  and  degradation 
mechanisms  of  the  individual  components  warrant  further  study. 
In  particular,  the  non-linear,  non-elastic  constitutive  behavior  of 
the  membranes  under  hygro-thermal  cycling  and  how  this  affects 
the  lifetime  of  the  fuel  cell  is  not  understood. 

In  a  previous  study,  we  investigated  the  stresses  in  the  PEM 
under  hygro-thermal  loading  (from  the  stress  free  condition  at 
assembly  to  operating  temperature  and  humidity)  [13,14].  How¬ 
ever,  in  that  study,  the  polymer  was  assumed  to  be  linear-elastic. 
Recent  experimental  investigations  in  our  laboratory  have  pro¬ 
vided  us  with  more  realistic  properties  of  PFSA-membranes 
[15].  The  experiments  show  that,  for  relatively  small  strains,  the 
mechanical  response  can  be  approximated  with  a  linear-elastic, 
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perfectly  plastic  constituent  behavior.  Thus,  linear-elastic,  per¬ 
fectly  plastic  material  response  with  temperature  and  humidity 
dependent  properties  will  be  incorporated  in  the  numerical  simu¬ 
lations.  Moreover,  we  will  investigate  the  response  of  the  mem¬ 
brane  during  one  cycle,  simulating  a  simplified  fuel  cell  duty 
cycle.  After  imposing  clamping  conditions  (described  below), 
this  cycle  consists  of  loading  (from  ambient  to  operating  condi¬ 
tion),  followed  by  unloading  (back  to  ambient,  corresponding  to 
shut-off).  We  will  see  that  significant  residual  stresses  develop 
during  unloading,  introducing  in-plane  tensile  stresses. 

In  the  following,  we  will  outline  the  assumptions  for  the 
material  behavior  followed  by  a  description  of  the  finite  element 
models  and  last  discuss  the  results  from  the  simulations. 


2.  Theory  for  isotropic  elasto-plasticity 


Here,  we  outline  the  approach  used  to  incorporate  hygro- 
thermal  effects  into  a  finite  element  simulation,  including 
isotropic  plasticity.  This  is  an  extension  of  previous  work 
[13,14],  where  only  the  elastic  behavior  of  the  membrane  was 
considered.  An  uncoupled  theory  is  assumed,  for  which  the  addi¬ 
tional  temperature  changes  brought  by  the  strain  are  neglected. 
We  assume  that  the  total  strain  tensor,  £y,  is  the  sum 


—  £?.  -f  £?!  -f  gT.  +  £?. 

IJ  '  ll  '  lJ  '  l]’’ 


where  £•)  is  the  elastic  strain  component,  sfj  the  plastic  strain 

component  and  sfj  and  sfj  are  the  temperature  and  swelling 
induced  strains,  respectively.  Let  7o  be  the  reference  temperature 
and  T  is  the  current  temperature.  The  thermal  strains  resulting 
from  a  change  in  temperature  of  an  unconstrained  isotropic  vol¬ 
ume  are  given  by 


sfj  =  a(T  -  T0)8jj  with  8 ij  = 


1  if  i  =  j 

0  if  i  ^  j  ’ 


where  a  is  the  linear  coefficient  of  thermal  expansion,  hav¬ 
ing  units  of  K_1,  and  8y  is  the  Kronecker  delta.  Similarly,  the 
swelling  strains  caused  by  moisture  uptake  are  given  by 
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where  C  is  the  relative  humidity  and  Co  is  the  reference  value 
for  C.  Moreover,  H)  is  the  temperature  and  humidity 

dependent  swelling-expansion  coefficient  tensor  due  to  mois¬ 
ture  absorption.  For  an  isotropic  material,  Pij(T,  H)  =  f$(T,  H)8y. 
The  swelling-expansion  coefficient  is  usually  defined  as  the  rel¬ 
ative  change  in  length  per  1%  water  weight  gain.  In  the  current 
application,  the  relative  humidity  rather  than  the  water  gain  is 
assumed  to  be  known,  thus  here  the  swelling-expansion  coeffi¬ 
cient  is  defined  as  the  relative  change  in  length  per  1  %  relative 
humidity  change.  Consequently,  the  units  of  are  1  %RH-1. 

Assuming  linear  response  within  the  elastic  region,  the 
isotropic  Hooke’s  law  is  used  to  determine  the  stress  tensor, 


E 


where  s\k  =  £n  +  £22  +  £33,  ^  Young’s  modulus  and  v  is 
Poisson’s  ratio.  Generalized  plane  strain  is  assumed  in  the  sim¬ 
ulations: 

£13  =  £31  =  £23  =  £32  =  0  and  £33  =  constant.  (5) 

For  the  inelastic  response,  incompressible  plastic  deforma¬ 
tion  is  assumed,  where  rate-independent  plastic  flow,  according 
to  the  von  Mises  yield  function  is  used  (J2-flow  theory)  [16]: 


fieri  j)  =  Y  -SySy  -  CTO,  (6) 

where  a y  is  the  components  of  the  (true)  stress  tensor  and  c>o 
is  the  yield  strength.  For  simplicity,  we  will  assume  that  the 
material  is  perfectly  plastic,  thus  the  yield  strength  does  not 
depend  on  the  plastic  strain  but  only  temperature  and  humidity. 
Therefore,  the  yield  strength  can  be  rewritten  as 

cr0  =  0-0(7: ;  H).  (7) 

Furthermore,  in  Eq.  (6)  Sy  are  the  components  of  the  deviatoric 
stress  tensor  is  defined  by 

1 

$ij  —  &ij  —&kk$ij-  (8) 

According  to  the  von  Mises  yield  criterion,  yield  occurs  when 

/K)  =  0.  (9) 

For/(<T/y)  <  0,  the  material  deforms  elastically. 

The  Mises  flow  theory  predicts  that  the  plastic  strain  incre¬ 
ment  tensor  must  be  proportional  to  the  deviatoric  stress  tensor, 
thus 

dsfj  =  SijdX,  (10) 

where  dX  is  a  scalar  proportionality  factor. 

3.  Material  properties 

The  material  properties  for  the  graphite  plates  are  set  to  those 
of  commercial  graphite  and  for  the  carbon  paper  from  TORAY® 
TGP-H-030  [18].  It  is  assumed  that  these  materials  have  linear- 
elastic  behavior  and  do  not  swell  in  response  to  moisture. 

Linear-elastic,  perfectly  plastic  constitutive  behavior  with 
temperature  and  humidity  dependent  material  properties  is 
assumed  for  the  membrane.  Young’s  modulus  and  yield  strength 
of  the  membrane  are  defined  for  four  temperature  and  relative 
humidity  values  based  on  our  experimental  data  for  Nation®  1 1 21 
[15].  The  tensile  tests  of  the  Nation®  films  were  conducted  in  an 
environmental  chamber  within  a  temperature  range  from  25  to 
85  °C,  and  a  relative  humidity  range  from  30  to  90%RH  to  obtain 
true  stress-true  strain  curves  at  each  temperature-humidity 
point  [15],  then  Young’s  modulus  and  yield  strength  are  cal¬ 
culated  at  these  points.  Even  though  a  slight  anisotropy  in  the 
material  properties  was  observed  experimentally  [15],  in  this 
study  we  assume  for  simplicity  that  the  material  properties  are 


(1  +  v)(l  -  2v) 


\v8ijslk  +  (1  -  2 v)sfj], 


(4) 


1  Nafion®  is  a  registered  trademark  of  E.I.  DuPont  De  Nemours  &  Co. 
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Table  1 

Material  properties  used  in  finite  element  analysis  [17,18] 


989 


Component  (material) 

k  (W  m-1  K-1) 

P  (kgm  3) 

E  (MPa) 

V 

a  (10-6  K-1) 

f  (mm mm  1  %RH  1 ) 

Cp  (Jkg-1  K-1) 

Bipolar  plates  (graphite) 

95 

1800 

10000 

0.25 

5 

Neglected 

750 

GDE  (carbon  paper) 

0.3 

400 

10000 

0.25 

-0.8 

Neglected 

500 

Membrane  (Nafion®  112) 

0.259 

2000 

Table  2 

0.25 

123 

f(T,  H)  (from  Fig.  1) 

1050 

isotropic.  Additional  physical  properties  for  the  membrane  were 
adopted  from  the  DuPont  Nation®  PFSA  Membranes  N-112 
Data  Sheet  [18].  The  material  properties  for  the  membrane  and 
other  components  used  in  this  study  are  shown  in  Tables  1-3  and 
Fig.  1. 


Table  2 

Young’s  modulus  at  various  temperatures  and  humidities  for  Nation®  112  [15] 


Young’s  modulusa  (MPa) 

Relative  humidity  (%) 

30 

50 

70 

90 

T=  25  °C 

197 

192 

132 

121 

U 

o 

in 

II 

161 

137 

103 

70 

T=  65  °C 

148 

117 

92 

63 

U 

o 

in 

oo 

II 

h-H 

121 

85 

59 

46 

a  Linear  interpolation  between  given  data  points  is  used  during  the  loading 
and  unloading. 


Table  3 

Yield  strength  at  various  temperatures  and  humidities  for  Nation®  112  [15] 


Yield  stressa  (MPa) 

Relative  humidity  (%) 

30 

50 

70 

90 

T=  25  °C 

6.60 

6.14 

5.59 

4.14 

U 

o 

in 

'n1‘ 

ll 

h-i 

6.51 

5.21 

4.58 

3.44 

r=65°C 

5.65 

5.00 

4.16 

3.07 

u 

o 

in 

oo 

II 

4.20 

3.32 

2.97 

2.29 

a  Linear  interpolation  between  given  data  points  is  used  during  the  loading 
and  unloading. 


Fig.  1 .  Experimental  data  for  dimensional  change  due  to  the  swelling  expansion 
and  the  polynomial  curve  fit  to  the  experimental  data  as  a  function  of  humidity 
and  temperature  plotted  for  four  constant  temperatures  for  Nation®  112  [15]. 


4.  Finite  element  model 

4.1.  Assumptions 

A  two-dimensional  finite  element  model  was  developed, 
assuming  generalized  plane  strain  conditions,  including  the  fol¬ 
lowing  simplifying  assumptions: 

(1)  Simplified  temperature  and  humidity  profiles  with  no  heat 
generation  were  assumed. 

(2)  The  electrodes  are  assumed  to  be  integrated  into  the  gas 
diffusion  layer  (GDL)  to  form  a  GDE  instead  of  considering 
this  as  a  separate  layer. 

(3)  For  the  GDE  and  the  bipolar  plate,  the  deformation  is 
assumed  elastic,  whereas  the  membrane  is  allowed  to 
deform  plastically. 

(4)  All  material  properties  are  assumed  to  be  isotropic  including 
the  expansions  due  to  the  thermal  load  and  moisture  uptake. 

(5)  One  hygro-thermal  load  cycle  (roughly  corresponding  to 
a  fuel  cell  duty  cycle)  is  applied  which  is  a  simultaneous 
loading  due  to  the  moisture  uptake  and  thermal  expansion 
followed  by  an  unloading  to  the  initial  state. 

The  model  used  in  the  current  investigation  is  an  extension 
of  a  model  developed  in  previous  studies  [13,14].  In  the  cur¬ 
rent  work,  we  incorporate  temperature  and  humidity  dependent 
properties  of  the  PEM,  as  well  as  including  plasticity  and  inves¬ 
tigating  the  response  to  one  simplified  duty  cycle. 

4.2.  Geometry 

A  schematic  drawing  of  a  PEMFC  unit  cell  used  in  the  FE 
model  is  shown  in  Fig.  2. 

Two  geometries  with  respect  to  the  gas  channels  are  modeled. 
The  gas  channels  in  the  graphite  plates  are  used  for  transport 
of  the  hydrogen  to  the  anode  side  of  the  membrane  electrode 
assembly  (MEA)  and  air  to  the  cathode,  while  water  vapor 
is  carried  away.  In  most  cases,  the  channels  on  either  side  of 
the  MEA  are  aligned  with  each  other.  However,  in  other  cases 
they  are  not  aligned,  for  example  due  to  manufacturing  flaws. 
Thus,  two  graphite  plate  alignments  are  selected  to  study  the 
effects  of  the  geometry  of  the  gas  flow  channels:  (1)  aligned 
and  (2)  alternating.  In  both  cases,  the  thickness  of  plates  is 
1 1  mm,  and  the  depth  and  the  width  of  grooves  in  the  plate 
are  1  mm. 

Furthermore,  two  clamping  methods  are  investigated:  (1) 
fixed  force,  corresponding  to  a  case  where  the  fuel  cell  stack  is 
equipped  with  springs  to  control  the  claming  force  and  (2 )  fixed 
displacement  where  the  graphite  plates  are  clamped  preventing 
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PEM  Fuel  Cell  Stack 


Fig.  2.  The  two  geometries  modeled,  where  the  flow  field  has  either  (A)  aligned  or  (B)  alternating  gas  channels.  The  mechanical  boundary  conditions  are  noted  in 
the  figure.  Two  additional  boundary  conditions  (not  shown)  for  case  (A)  and  (B)  were  imposed  by  subjecting  the  stack  to  either  constant  displacement  or  constant 
force  clamping  at  the  top. 


any  overall  deformation  in  the  thickness  direction.  In  the  former 
case,  a  constant  pressure  (1  MPa)  is  applied  on  the  surface  of 
upper  graphite  plate.  In  the  latter  case,  the  applied  displacement 
is  calculated  by  applying  a  1  MPa  load  to  an  otherwise  unloaded 
model,  and  the  resulting  displacement  at  the  upper  boundary  is 
fixed  throughout  the  analysis. 

The  membrane  thickness  is  50  jim  (corresponding  to  the 
thickness  of  Nation®  112).  As  mentioned  previously,  the  cata¬ 
lyst  layer  is  included  in  the  gas  diffusion  layer,  forming  a  gas 
diffusion  electrode  (GDE)  with  a  thickness  of  100  pan. 

The  commercial  software,  ABAQUS  6.4  [19]  is  used  to 
conduct  the  simulations.  To  simulate  a  unit  cell,  continuous 
boundary  conditions  are  applied  on  the  model  (Fig.  2).  These 
are  manifested  as  uz  =  0  applied  on  the  bottom,  and  ux  =  0  on 
the  right  side  (Fig.  2).  A  linear  constraint  A\u =  0 
is  applied  along  the  left  edges  of  the  unit  cell.  Here,  A\  -  1  and 
A2  =  —  1 ,  N\  is  the  node  at  the  bottom  left  and  N2  are  the  remain¬ 
ing  nodes  on  the  left  side.  This  condition  constrains  the  right 


hand  side  to  displace  uniformly.  Four-node  generalized  plane 
strain  temperature-displacement  coupled  elements  (CPEG4T) 
are  used  in  the  FE  analysis. 

4.3.  Simulating  eigenstrain:  swelling  and  thermal 
expansion 

The  total  strain  in  the  membrane  due  to  the  moisture  and  tem¬ 
perature  change  is  calculated  from  the  coefficients  of  swelling 
and  thermal  expansion.  To  achieve  this  in  the  simulations,  the 
user  subroutine  UEXPAN  [19]  is  used.  A  constant  thermal 
expansion  coefficient,  a  (Table  1),  and  a  swelling-expansion 
coefficient,  /3(T,  //),  expressed  as  a  polynomial  function  of 
humidity  and  temperature  are  used  in  this  subroutine  to  define 
the  expansion  of  the  membrane  from  initial  state  (35%RH  and 
20  °C)  to  the  hydrated  state  (100%RH  and  85  °C).  UEXPAN 
calculates  an  eigenstrain  based  on  the  state  of  the  defined  field 
variables  humidity  and  temperature,  relative  to  the  initial  state 
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Table  4 

Constants  of  the  swelling  strain  polynomial  (Cy;)  as  defined  in  Eq.  (11) 


Temperature  (°C) 


991 


i=  1 

i  =  2 

i  =  3 

i  =  4 

7  =  1 

3.0816E— 12 

— 4.4727E— 10 

2.2696E— 08 

— 7.5602E— 08 

7  =  2 

— 5.3516E— 10 

7.6152E— 08 

— 3.7261E— 06 

2.3356E— 05 

7  =  3 

3.6359E— 08 

— 5.2865E— 06 

2.6400E— 04 

— 2.1601E— 03 

7  =  4 

— 6.9232E— 07 

1.0213E— 04 

— 5.1794E— 03 

4.5824E— 02 

(recall  that  the  materials  are  assumed  to  have  isotropic  thermal 
and  swelling-expansion  behavior  that  is  uncoupled). 

The  coefficient  of  swelling  expansion,  f3(T ,  //),  is  derived 
from  the  swelling  strain  polynomial  function  which  is  obtained 
by  curve  fitting  to  the  experimental  data  of  swelling  strains 
(Fig.  1).  The  swelling  strain  polynomial  is  defined  as 

4 

ss(T,  H)  =  Cjj  T4~‘  RH4-;  ,  (1 1) 

ij=  1 

where  T  and  H  are  the  temperature  (°C)  and  relative  humidity 
(%),  respectively.  The  constants  of  the  polynomial  Cp  are  given 
in  Table  4. 

4.4.  Simulation  of  a  fuel  cell  duty  cycle 

A  fuel  cell  duty  cycle  is  simulated  by  first  defining  “initial 
conditions,”  followed  by  conducting  three  (ABAQUS)  steps: 

The  initial  conditions  corresponding  to  zero  stress-state  are 
defined;  all  components  of  the  cell  stack  are  set  to  room  tem¬ 
perature  20  °C,  and  relative  humidity  35%  (corresponding  to 
the  assembly  conditions).  Since  we  do  not  have  experimen¬ 
tal  data  for  temperatures  below  25  °C  and  relative  humidities 
below  35%,  the  material  properties  used  below  these  values  are 
assumed  to  be  constant. 

•  Step  T.  Impose  the  clamping  condition  (fixed  force  or  fixed 
displacement  as  described  above). 

•  Step  2:  Hygro-thermal  loading  is  applied  by  linearly  increas¬ 
ing  the  relative  humidity  of  the  membrane  up  to  100%,  and 
defining  an  estimated  temperature  profile  as  follows:  cath- 
ode/GDE  and  anode/GDE  interface  temperatures  are  at  86 
and  85  °C,  respectively,  and  the  bipolar  plate  is  assumed  to 
have  a  constant  temperature  at  its  mid-plane,  80  °C  which  is 
the  cell  operating  temperature.  These  values  are  reasonable 
for  a  fuel  cell  operating  at  80  °C  under  typical  conditions 
[15].  During  the  cycle,  the  internal  temperature  distribution 
profile  in  the  membrane  was  calculated  from  these  bound¬ 
ary  conditions.  The  resulting  temperature  profile  is  shown  in 
Fig.  3. 

•  Step  3:  Unloading  is  achieved  by  linearly  decreasing  the  rela¬ 
tive  humidity  and  temperature  to  the  initial  relative  humidity 
and  temperature  values,  but  leaving  the  clamping  conditions 
imposed. 

Thus,  this  simplified  duty  cycle  corresponds  to  taking  the  cell 
from  start-up  to  hydrated  operation  and  then  back  to  shut-down, 
assuming  the  cell  is  left  to  dryout  during  shut-down. 


80  85  86 


Fig.  3.  The  temperature  profile  for  the  unit  fuel  cell,  assuming  a  constant 
anode/GDE  temperature  of  85  °C  and  cathode/GDE  temperature  of  86  °C 
[13,14]. 

5.  Results  and  discussion 

5.1.  Stress  and  deformation  history  of  the  selected  points 

First,  we  will  investigate  the  stress  histories  at  two  selected 
points  of  the  membrane  for  each  case,  throughout  the  cycle. 
The  “left  end”  and  “right  end,”  as  marked  in  Fig.  2,  are 
tracked  to  elucidate  the  evolution  of  the  stresses  through  the 
loading-unloading  sequence.  The  effective  stresses  according 
to  von  Mises,  “Mises  stresses,”  given  by 


crM  =  ]J-sijsij  (12) 

are  a  useful  measurement  of  the  effective  stress.  When  <tm  >  <to> 
yielding  occurs,  according  to  Eq.  (6). 

The  Mises  stresses  that  develop  during  the  cycle  are  plot¬ 
ted  together  with  the  yield  strength  in  tension  and  compression 
to  determine  when  the  yield  strength  is  reached  and  the  plastic 
deformation  begins  (Fig.  4).  Yield  strength  for  each  temperature 
and  humidity  combination  reached,  was  determined  using  a  lin¬ 
ear  interpolation  method  between  the  values  of  the  yield  strength 
obtained  at  16  temperature-humidity  points  from  the  experi¬ 
ments  [15].  Although  the  Mises  stress  is  by  definition  always 
positive,  we  have  associated  a  sign  of  the  Mises  stress  with  the 
trend  of  the  major  stress  component  (in-plane  stress)  for  com¬ 
parison.  This  way  a  general  trend  can  be  determined  revealing 
overall  compressive  or  tensile  behavior.  The  history  plots  for  the 
in-plane  stresses  are  shown  in  Fig.  5.  From  these  results,  it  can 
be  seen  that  the  Mises  stress  strongly  depends  on  the  in-plane 
stress  component. 
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Fig.  4.  The  evolution  of  the  Mises  stress  during  one  simulated  fuel  cell  duty  cycle  plotted  with  the  yield  strength  of  the  membrane  for  the  selected  points  for  (A) 
aligned,  fixed  displacement;  (B)  aligned,  fixed  force;  (C)  alternating,  fixed  displacement  and  (D)  alternating,  fixed  force.  The  apparent  sign  of  the  Mises  stress  is 
obtained  by  following  the  trend  of  the  in-plane  stress.  The  cycle  is  characterized  with  hygro-thermal  loading  from  initial  state,  20  °C,  35%RH,  to  maximum  load 
85  °C,  100%RH,  and  for  unloading  from  maximum  load  back  to  initial  state. 


Initially  the  membrane  is  pre-stressed  due  to  the  applied 
clamping  as  explained  previously.  During  the  hygro-thermal 
loading,  the  in-plane  stress  and  the  Mises  stress  increase  in 
compression  to  eventually  reach  yielding  (Fig.  4).  After  initial 
yielding  occurs  the  internal  forces  are  redistributed  to  accom¬ 
modate  additional  loading.  Upon  unloading,  the  compressive 
stresses  decrease,  eventually  leading  to  tensile  stresses  (Fig.  4). 
Thus,  due  to  the  yielding  at  elevated  temperatures,  the  load  cycle 
imposes  a  residual,  tensile  stress  in  the  membrane. 

5.2.  Stress  distribution  along  the  membrane 

With  the  evolution  of  the  stresses  (for  selected  points) 
throughout  the  cycle  established,  we  will  next  investigate  the 
overall  stress  distribution  at  maximum  temperature  and  after 
unloading.  Fig.  6  shows  the  stress  distributions  in  the  membrane 
for  the  alternating — fixed  force  case  at  maximum  hygro-thermal 


load  and  after  unloading.  At  maximum  load  (85  °C,  100%RH), 
the  membrane  is  in  overall  yielding  (Fig.  6A).  As  discussed  in 
the  previous  sub-section,  the  yielding  results  in  a  redistribution 
of  the  internal  forces  and  stresses.  Consequently,  tensile  in-plane 
stresses  are  achieved  upon  unloading  (20  °C,  35%RH)  (Fig.  6B). 
We  note  that  if  yielding  did  not  occur  during  loading,  the  stresses 
after  unloading  would  be  identical  to  the  initial  stresses  (stresses 
due  to  the  clamping  -  fixed  force  or  fixed  displacement  -  only). 

Fig.  7  shows  the  stress  distribution  along  the  top  edge 
of  the  membrane,  after  loading  and  unloading  for  the 
alternating — fixed  displacement  case.  Considering  the  in-plane 
stress  and  the  Mises  stress  distribution  for  the  alternating  case 
(Figs.  8B  and  9B)  suggest  that  the  Mises  stress  is  dominated 
by  the  in-plane  stress  during  the  deformation.  The  out-of-plane 
stress  is  always  smaller  than  the  in-plane  stress,  although  they 
are  at  the  same  order  of  magnitude  in  some  limited  situations.  In 
the  confined  regions  of  the  membrane  (i.e.  midpoint  of  the  mem- 
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Fig.  5.  The  evolution  of  the  in-plane  stress  during  one  fuel  cell  duty  cycle  plotted  for  the  selected  points  for  (A)  aligned  and  (B)  alternating  case  for  each  clamping 
method.  The  cycle  is  characterized  with  hygro-thermal  loading  from  initial  state,  20  °C,  35%RH,  to  maximum  load  85  °C,  100%RH,  and  for  unloading  from  maximum 
load  back  to  initial  state. 


(A)  Mlses  Stress  (MPa) 


Stresses  at  Maximum  Hygro-Thermal  Load  ( 100%RH.  85°C  ) 


Residual  Stress  after  Unloading  ( 35%RH.  20°C ) 


(B)  In-Plane  Stress  (MPa) 


/  /  I 

Stresses  at  Maximum  Hygro-Thertnal  Load  ( 100%RH.  85°C ) 
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(C) Shear  Stress  (MPa)  r  - 


Stresses  at  Maximum  Hygro-Thermal  Load  ( 100%RH.  85°C  ) 
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Fig.  6.  The  contour  plots  of  the  membrane  for  alternating — fixed  displacement,  at  maximum  hygro-thermal  load  and  after  unloading  (“residual  stress”),  (A)  Mises 
stress,  (B)  in-plane  stress,  (C)  in-plane  shear  stress  and  (D)  out-of-plane  stress. 
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Fig.  7.  Shear  stress  and  the  out-of-plane  stress  at  maximum  hygro-thermal  load 
and  after  unloading  for  alternating,  fixed  displacement  as  a  function  of  location 
along  the  membrane. 


Fig.  9.  Mises  stress  at  maximum  hygro-thermal  load  (85  °C,  100%RH)  and  after 
unloading  (20  °C,  35%RH)  for  (A)  aligned  and  (B)  alternating  cases  and,  two 
different  clamping  methods  as  a  function  of  location  along  the  membrane.  <t({ 
and  <7q  denote  the  yield  strengths  at  20  °C,  35%RH  and  at  85  °C,  100%RH, 
respectively. 

brane  in  Fig.  7),  the  out-of-plane  stress  becomes  large  enough 
to  have  the  same  contribution  as  the  in-plane  stress.  This  is  par¬ 
ticularly  true  in  the  case  of  fixed  displacements.  Shear  stresses 
are  small  and  may  be  neglected  for  all  other  cases.  The  in-plane 
and  Mises  stress  distribution  for  the  aligned  case  are  shown  in 
Figs.  8 A  and  9 A  for  comparison. 

The  plastically  deformed  regions  can  be  visualized  from  the 
Mises  stress  distribution  (Fig.  9A  and  B).  The  regions  of  the 
membrane  that  deform  plastically  after  hygro-thermal  loading 
depend  on  the  clamping  method  and  the  alignment  of  the  model. 
For  the  aligned-fixed  displacement  case,  the  right  part  of  the 
membrane  (Fig.  9 A)  remained  elastic  during  the  deformation. 
In  this  case,  the  right  part  of  the  membrane  is  constrained  from  all 
sides  which  leads  to  large  hydrostatic  compression  component.2 


Fig.  8.  In-plane  stress  at  maximum  hygro-thermal  load  (85  °C,  100%RH)  and  - 

after  unloading  (20  °C,  35%RH)  for  (A)  aligned  and  (B)  alternating  cases,  for  2  According  to  Mises  yield  criterion,  Eq.  (9),  hydrostatic  pressure  does  not 
two  clamping  methods  as  a  function  of  location  along  the  membrane.  influence  yielding. 
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Fig.  10.  Maximum  values  of  the  stress  components;  in-plane,  out-of-plane  and 
shear,  respectively,  induced  in  the  membrane  by  the  hygro-thermal  loading  are 
shown  for  each  alignment  and  clamping  method.  Letters  on  the  columns  indicate 
the  location  of  the  stress  in  the  membrane,  where  L  corresponds  to  the  left  end; 
M  to  the  midpoint,  and  R  to  the  right  end. 

For  the  aligned  case,  the  right  end  of  the  membrane  is  con¬ 
strained  between  the  GDE  and  bipolar  plates,  and  the  left  end  of 
the  membrane  can  expand  freely  due  to  the  thermal  and  swelling 
strains,  which  eventually  leads  to  different  deformation  behavior 
and  stress  values.  For  the  alternating  case,  due  to  the  symmetry 
with  respect  to  midpoint  of  the  membrane,  the  left  and  right  ends 
have  the  same  stresses  and  the  middle  point  has  the  maximum 
stresses. 

5.3.  Maximum  stresses  in  the  membrane 

Finally,  we  will  consider  the  magnitude  of  the  maximum 
stresses.  The  maximum  stresses  generated  in  the  membrane  dur¬ 
ing  the  hygro-thermal  cycle  are  summarized  in  Fig.  10.  The 
location  of  the  maximum  stress  is  indicated  and  classified  as 
belonging  to  the  left  side  (L),  the  middle  (M),  or  the  right 
side  (R)  of  the  membrane.  The  locations  are  determined  from 
Figs.  7  and  8  according  to  the  following: 

(i)  As  indicated  in  Fig.  10,  the  in-plane  stresses  are  always 
higher  than  the  out-of-plane  stresses.  However,  for  fixed 
displacement,  the  out-of-plane  stresses  are  close  to  the  in¬ 
plane  stresses  and  their  maxima  are  at  the  same  points:  at 
the  right  end  for  the  aligned  case,  and  at  the  midpoint  for 
the  alternating  case. 

(ii)  For  the  fixed  force  case,  the  maximum  in-plane  stresses 
are  tensile,  whereas  the  out-of-plane  stresses  are  smaller 
and  compressive.  Since  the  fixed  force  clamping  allows 
expansion  of  the  unit  cell,  the  membrane  is  not  con¬ 
strained  along  the  grooves  and  tensile  stresses  develop 
upon  unloading.  Therefore  the  maximum  in-plane  stresses 
occur  at  the  left  end  for  the  aligned  case,  and  at  the  left 
and  right  end  for  the  alternating  case  (Fig.  8).  The  max¬ 
imum  of  the  out-of-plane  stresses  are  at  the  right  end  for 
the  aligned  case,  and  at  the  midpoint  for  the  alternating 
case. 


Fig.  11.  Maximum  in-plane  residual  stresses  obtained  at  the  three  points  of  the 
membrane;  the  left  end,  the  midpoint  and  the  right  end,  shown  for  each  alignment 
and  clamping  method. 

(iii)  For  all  cases,  the  shear  stresses  are  small,  with  the  alternat¬ 
ing  case  causing  the  higher  shear  stresses  in  the  membrane. 
The  maximum  shear  stresses  are  at  the  midpoint  of  the 
membrane  for  all  cases  (Fig.  7). 

The  residual,  in-plane  stresses  found  after  unloading  are  plot¬ 
ted  at  the  three  points  (Left,  Middle,  Right)  of  the  membrane 
which  correspond  the  maximum  and  the  minimum  values  of  the 
stresses.  As  seen  from  Fig.  1 1,  the  residual  in-plane  stresses  are 
tensile,  except  for  the  right  end  of  the  membrane  in  the  aligned- 
fixed  displacement  case  since  no  plastic  deformation  occurs,  as 
discussed  in  Section  5.2.  Also,  for  all  cases  the  maximum  resid¬ 
ual  stresses  are  similar. 

6.  Conclusions 

The  mechanical  response  of  a  PEM  subjected  to  a  single 
hygro-thermal  loading  cycle,  simulating  a  simplified  single  fuel 
cell  duty  cycle,  was  investigated.  A  linear,  uncoupled  tempera¬ 
ture  and  humidity  profile,  assuming  steady-state  conditions,  was 
used  for  the  loading-unloading.  Linear-elastic,  perfectly  plas¬ 
tic  material  response  with  temperature  and  humidity  dependent 
material  properties  was  used  to  study  the  plastic  deformation 
behavior  of  the  membrane  during  the  cycle.  Independent  of  the 
clamping  method  and  alignment,  plastic  deformation  always 
occurs.  However,  the  residual  stresses  in  the  membrane  directly 
correlate  with  alignment  and  clamping  method. 

For  “fixed  displacement stresses  are  generally  higher  than 
those  assuming  “fixed  force."  The  shear  stresses  are  small  for 
all  cases,  although  alternating  gas  channel  alignment  produces 
higher  shear  stresses  (less  than  15%  of  the  in-plane  stress  at 
maximum).  The  out-of-plane  stress  is  small  compared  to  the 
in-plane  stress,  especially  for  the  fixed  displacement  case.  In 
general,  the  in-plane  stress  is  the  largest  stress  component  during 
the  loading  sequence,  thus  it  controls  the  yielding  behavior. 

Hygro-thermal  loading  causes  compressive  stresses  in  the 
membrane  due  to  the  constraints,  which  prevent  expansion  due  to 
the  swelling  of  the  membrane.  During  unloading,  the  membrane 
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shrinks  and  due  to  redistribution  of  the  stresses  caused  by  the 
plastic  deformation,  tensile  residual  stresses  may  develop  dur¬ 
ing  cooling/drying.  In  many  cases,  the  residual  stresses  become 
large  enough  to  cause  tensile  yielding.  We  believe  that  these 
high  in-plane  residual  stresses  are  a  significant  contributor  to 
the  mechanical  failures  observed  in  the  membranes  since  they 
may  cause  the  propagation  of  through-the-thickness  flaws  or 
tearing  in  the  membrane. 

In  conclusion,  we  have  shown  that  use  of  temperature  and 
humidity  dependent  material  properties  in  the  mechanical  mod¬ 
eling  of  the  membranes  leads  to  critical  residual  stresses  after  one 
fuel  cell  duty  cycle.  The  results  suggest  that  the  in-plane  residual 
tensile  stresses  developed  upon  unloading,  may  lead  to  the  fail¬ 
ure  of  the  membranes  due  to  the  mechanical  fatigue.  Therefore, 
in  order  to  acquire  a  complete  understanding  of  these  damage 
mechanisms  in  the  membranes,  mechanical  response  under  con¬ 
tinuous  hygro-thermal  cycles  should  be  studied  under  realistic 
cell  operating  conditions. 
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